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Abstract 

Nuclear processes in solid environment are investigated. It is shown that if a slow, quasi-free 
heavy particle of positive charge interacts with a "free" electron of a metallic host, it can obtain 
such a great magnitude of momentum in its intermediate state that the probability of its nuclear 
reaction with an other positively charged, slow, heavy particle can significantly increase. It is also 
shown that if a quasi-free heavy particle of positive charge of intermediately low energy interacts 
with a heavy particle of positive charge of the solid host, it can obtain much greater momentum 
relative to the former case in the intermediate state and consequently, the probability of a nuclear 
reaction with a positively charged, heavy particle can even more increase. This mechanism opens 
the door to a great variety of nuclear processes which up till know are thought to have negligible 
rate at low energies. Low energy nuclear reactions allowed by the Coulomb assistance of heavy 
charged particles is partly overviewed. Nuclear pd and dd reactions are investigated numerically. 
It was found that the leading channel in all the discussed charged particle assisted dd reactions is 
the electron assisted d + d — ?■ "^He process. 

PACS numbers: 25.70.Jj, 25.45.-z, 25.40.-h 

Keywords: fusion and fusion-fission reactions, ^i/-induced nuclear reactions, nucleon induced reactions 



* retired from Budapest University of Technology and Economics, Institute of Physics, 
e-mail: kalmanpeter3@gmail.com 



1 



I. INTRODUCTION 



It is a standard of nuclear physics that heavy, charged particles j and k of like positive 
charge of charge numbers zj and z^. need considerable amount of relative kinetic energy E 
determined by the height of the Coulomb barrier in order to let the probability of a nuclear 
interaction have significant value. Mathematically it appears in the energy dependence of 
the cross section (a) of the charged-particle induced reactions as 

aiE) = S {E) exp [-2Tii^,k {E)] /E, (1) 

where S (E) is the astrophysical factor, which can be written as S{E) = S{0) + SiE + S2E'^ , 
The Sommerfeld parameter 

Vjk{k) = ZjZkOf^^, (2) 

where k = |kj — k/d is the magnitude of the relative wave vector k = kj— k^ of the interacting 
particles of wave vectors kj and k^, (k ~ \/E). The reduced mass of particles j and k of 
rest masses rrij and rrik 

fijk = rrijTnk/ {rrij + rrik) , (3) 

h is the reduced Planck constant, e is the elementary charge and a/ is the fine structure 
constant. (It can be shown that in the case of slow relative motion the exponential function 



in ([T]) is the same as the Gamow factor |2|, that hinders nuclear reactions between particles of 
like electric charge.) The energy dependence of the cross section ([T]) can be derived applying 
the Coulomb solution 

^(r) = e^''-7(k,r)/v^, (4) 

which is the wave function of a free particle of charge number zj in a repulsive Coulomb field 
of charge number z^, where V denotes the volume of normalization, and r is the relative 
coordinate of the two particles. Here 

/(k, r) = e-^^'^^/^r^^ ^ tr]jk)iFi{-tr],k, l;i[kr-k- r]), (5) 



where iFi is the confluent hypergeometric function and T is the Gamma function 

It is the consequence of energy dependence ([T]) of the cross section that to this day it 
is a commonplace that the rate of any nuclear reaction between heavy, charged particles 
of positive charge is unobservable at low energies. The aim of this paper is to show that 



in a solid (particularly in a metal), contrary to the former assumption, there are nuclear 
processes that can have observable rate at low energies. 

II. PRELIMINARY CONSIDERATIONS 

In the low-energy range {kR ^ 1, where R is the radius of a nucleon) and for |r| < i? 
the long wavelength approximation 

|y.(r)| = 1^(0)1 = /,,(A:)/v^ (6) 

is valid, where 

is the Coulomb factor. We introduce the notation 

FAk) = ffkik) (8) 

with which the cross section (and the rate) of a first order process is proportional. If k is 
small then the long wavelength approximation produces the form ([T]) of cr well. Thus the 
fact that the rate of any nuclear reaction between heavy, charged particles of positive charge 
is unobservable at low energies is the consequence of Fjk{k) being small. 

In solids, however, where free electrons are present nuclei also Coulomb interact with 
electrons. It means that the Hamiltonian governing the state of the nuclei must also contain 
the interaction Hamiltonian with the electrons. If one wants to describe the interaction 
(scattering) of a nucleus with an other one that takes into consideration the interaction 
with the electrons the lower order process appears in the perturbation calculation can be 
seen in FIG. 1. The basic idea can be demonstrated with the aid of FIG. 1(a), in which a 
Coulomb scattering is followed by a capture process governed by strong interaction. When 
calculating the transition probabihty (and the rate) of such a second order process the 
following statements are valid. Energy and momentum (wave number vector) are conserved, 
i.e. Ei = Ef, kj = kj, where Ei and Ef are the total energies, and kj and kj are the 
total wave number vectors in the initial and final states, respectively. However energy-wave 
number vector (momentum) conservation may be violated in the "intermediate" state. In 
the cases investigated the initial particles (particles 1, 2 and 3) are slow and the sum of their 
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initial kinetic energies Ei and the sum of their wave number vectors kj can be neglected, i.e. 
Ei = and kj = can be supposed. 

It is thought that particle 1 is an electron and particles 2 and 3 are heavy, and of positive 
charge. The nuclear reaction 2 + 3 — 4 has reaction energy A. This energy is shared 
between the outgoing particles 1 and 4. Thus particle 1 obtains energy and wave number 
vector of nuclear order of magnitude. Since the Coulomb interaction in the case of free 
particles conserves wave number vector (momentum), and since the initial wave number 
vector of particles 1 and 2 can be neglected in wave number conservation, particle 2' gets 
a wave number vector k2/ opposite to the final wave number vector ki/ of particle 1', i.e. 
k2' = — ki/. Moreover if one calculates the Coulomb matrix element using plane waves for 
the freeparticles then the matrix element must be corrected with the so called Sommerfeld 
factor jj| 

_ /l2(|k2-ki|) 

where /12-s are Coulomb factors [see ([7])] for particles 1 and 2 of electric charge numbers 
Z2 = 1 and Zi = —1 since particle 1 is an electron. If particle 2 is heavy and slow, and particle 
1 is an electron then |k2 — ki| = /ci, that is the magnitude of the initial wave number vector 
of particle 1, furthermore |k2' — ki/| = 2/c2' since k2' = — ki/. Thus the cross section and the 



rate are proportional to 



Gs{h,2k,,)=gl = p^ (10) 

It can be shown that ^12(2^2') = 1 and therefore Gs = Fuih) = 23.18 x (^E^^^^ {eV)^ 
where Ei is the energy of the initial free electron. When calculating the matrix element 
of the strong interaction potential between particles 2 and 3 we use (EI). Consequently 
the second order rate is proportional to GsE23 which is mainly determined by -^23(^2') — 
exp [— 2717723 {k2')]- Now k2> = |ki/| has a nuclear order of magnitude. 

The first order rate is proportional to F2z{k) ~ exp [—2717723 (fc)], where k is the wave 
number of the slow, initial particle 2 (particle 3 is supposed to beat rest). Since k <^ k2', 
and therefore 7723 (fc) ^ V23{k2'), and consequently exp [—2717723 (/c)] ^ exp [—2717723 (A;2')], 
the rate of the second order process is much higher than that of the first order process. As 
a result, although the rate of a second order process is usually much less than the rate of a 
first order process, in this case the exponential increment is so huge that it can dwarf the 
rate of the first order process. 
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As a numerical example we consider the electron assisted d + d ^ process with slow 
deuterons. In this case, one of the slow deuterons (as particle 2) can enter into Coulomb 
interaction with a quasi-free, slow electron of the solid before the nuclear reaction (see FIG. 
1(a)). The states of the free deuteron and the free electron can be described by plane waves, 
therefore the Coulomb interaction preserves the wave number vector (momentum). If in 
the second order process the Coulomb interaction is followed by strong interaction, which 
induces a nuclear capture process, then the energy A of the nuclear reaction is divided 
between the electron and the heavy nuclear product. Since the rest mass rriN oi the nuclear 
product is much larger than the rest mass mg of the electron, the electron will take almost all 
the total nuclear reaction energy A away and the magnitude of its wave number vector A^/ = 
\/A2 + 2meC^A/ (he) ~ A/ (he) (if A 3> nieC^)- If initially (before the Coulomb interaction) 
the electron and the deuteron move slowly and the magnitudes of their wave number vectors 
are much smaller than A/ (He), then the initial wave number vectors can be neglected in 
the wave number vector (momentum) conservation and consequently, in the intermediate 
state the deuteron gets a wave number vector of magnitude kr and of direction opposite 
to the wave number vector of the electron. Thus the deuteron will have a (virtual) wave 
number vector in the intermediate state that is large enough to make it able to overcome 
the Coulomb barrier by tunneling and to take part in a nuclear process. If A = 23.84 
[Mey] , which is the reaction energy of the d + d ^ ^-^^e reaction, then the deuteron will 
have a (virtual) wave number k2' = A/ (he) in the intermediate state (in state 2'). The 
corresponding value F23 = 0.356. It must be compared e.g. to the extremely small value 
F23 (1 eV) — 1.1 X 10~^^^ that is characteristic of the first order process. 

III. ELECTRON ASSISTED NUCLEAR PROCESSES 

The change of state of heavy charged particles induced by solid state environment is 
modelled in the following way. Let us take two independent systems A and B, where A is a 
solid and B is an ensemble of free, heavy charged particles (e.g. a free deuteron or proton 
gas) with the corresponding Hamiltonians Ha and Hb- It is supposed that their eigenvalue 

problems arc solved, and the complete set of the eigenvectors of the two independent systems 
are known. Let us extend the state vectors of systems A and B to those nuclear bound states, 
which are initially empty, corresponding to the assumption that at the beginning the two 
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systems do not interact. The interaction between them to be switched on adiabatically is 
described by the interaction Hamiltonian Vab = (xab) + V^^ (xab), where V'-^'^ and V^* 
stand for the Coulomb and the strong interaction potentials, respectively, and the suffixes A 
and B in their argument symbolize that one party of the interaction comes from system A 
and the other from system B. (Similar model is used by jsl introducing the reduced density 
operator.) In the process investigated, ffist a heavy, charged particle of system B takes 
part in a Coulomb scattering with any charged particle of system A and it is followed by a 
strong interaction with some nucleus of system A that leads to their final bound states. The 
graphs of the process can be seen in FIG. 1. The process of FIG. 1(a) is a nuclear capture 
process and the process of FIG. 1(b) is a nuclear reaction. (The processes, where the nuclear 
interaction is followed by the Coulomb interaction, can be neglected because of the situation 
discussed in the introduction.) Particles 1 and 3, belonging to system A are: electrons (e.g. 
free electrons in the case of a metal) and localized heavy, charged particles (bound, localized 
p, d and other nuclei) as nuclear targets. Particle 2 belongs to system 5, that is a charged, 
heavy particle (e.g. proton (p) or deuteron (rf)), and that is supposed to move freely in a 
solid (e.g. in a metal). Since the aim of this paper is to show the fundamentals of the main 
effect, the problem, that there may be identical, indistinguishable particles in systems A 
and B is not considered, the simplest description is chosen, and the dynamic evolution of 
the number N2 of particles 2 of system B is not investigated. 

The transition probability per unit time ^W^/^'* j of the process can be written as 

< = |EKt^(^/-^) (11) 

^ f 

with 

^^?' = Ef^^*(k./ + K„„,). (12) 

where Ei, and Ej are the kinetic energies in the initial, intermediate and final states, 
respectively, A is the reaction energy, i.e. the difference between the rest energies of the 
initial and final states, V is the volume of normalization. is the matrix element of 

the Coulomb potential between the initial and intermediate states and V^^* is the matrix 
element of the potential of the strong interaction between intermediate and final states. 
type = a, +P,nP correspond to the reaction of FIG. 1(a) (a) and to the reactions of FIG. 
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(b) 



FIG. 1: The graphs of electron assisted nuclear reactions. The simple lines represent free (initial 
(1) and final (1')) electrons. The doubled lines represent free, heavy, charged initial (2) particles 
(such as p, d), their intermediate state (2'), target nuclei (3) and reaction products (4, 5, 6). The 
filled dot denotes Coulomb-interaction and the open circle denotes nuclear (strong) interaction. 
FIG. 1(a) is a capture process and FIG. 1(b) is a reaction. 

1(b) with both particles charged (+/3) and with one of them neutral (n/3), respectively. 

= Elf + El^type), (13) 



E, = Elf (ki^) + E, (k2^) , (14) 

where Eif is the kinetic energy and kij is the wave vector of particle 1 (electron) in the 
final state. E(^a) = -E'4 (k4), K(q,) = k4 in the case of process of FIG. 1(a) and or np) = 
E^ (ks) + Eq (ke), K(+/3 or np) = ks + kg in the case of process of FIG. 1(b) with Ej (kj) the 
kinetic energy of the j — th particle in the intermediate (particle 2) and final (particles 4 or 
5, 6) states. Since particle 1 is an electron, 

Elf = ^{hcfklf + mlc^ - rriec^ (15) 

with nieC^ denoting the rest energy of the electron. 
For the Coulomb potential we use its screened form 

V^\-k) = j ^^^exp(.q.x)dq (16) 
7 



with screening parameter A and coupling strength = afhc. For the strong interaction the 
interaction potential 

V'\.) = -2f'-^^^i^ (17) 

is applied, where the strong coupling strength = O.OS^c [6|] and 1/s is the range of 
the strong interaction. The calculation of the total rates of the electron assisted nuclear 
processes can be found in Appendix I. 

The result of the total rate of the leading, electron assisted {p or d) capture process 

Wi^i {a) = Ktot {a) {Gs{ku, 2hf))^^ B^!^hl^^^^uN, (18) 

with kif = ^. Here u denotes the deuteron (or proton) over metal number densities, N2 
is the number of initial particles 2 and for Gs see ffTOj) . Furthermore, we introduced the 
following notation 

32 

Ktot (a) = -^geKo («) , (19) 
where Qe is the number of the valence electron states corresponding to one unit cell and 

(a) = 192 (277)^ zlzl - a) {Hcf Rc. (20) 

Here R = 1.2 x lO^^'^ [cm] is the radius of a nucleon (corresponding to the single nucleon 
approach applied). Finally, hcorr,3 = — in the case of proton capture process and 
hcorr,3 = ^3 in the case of deuteron capture reactions (both are taken in the Weisskopf 
approximation), where A3 and Z3 are the mass and charge numbers of particle 3. 

Taking Z2 = 1, a/ = 1/137, f/ {he) = 0.08, c = 3 x 10^° [cm/s], Vc = d^/A (the volume 
of unit cell for Ni and Pd) with d = 3.89 x 10^^ [cm] {Pd lattice) one gets Ktot («) = 560 
[MeV^s-^] in the case of Pd. Furthermore {Gs)^^ = {Fuihi))^^ = 23.18 x (e~^^^ (eV)) , 

\ / av 

where Ei is the energy of the initial free electron and the average is made by means of Fermi- 
Dirac distribution. (A^2 and F23{kif = ■^) are dimensionless and A has to be substituted in 
[MeV] units in (fTSj).) 

The quantity F2z{kif = ■^) given by is the most important factor in Wt^t («)• The 
electron (particle 1) transfers large momentum (energy) to particle 2' through Coulomb 
scattering. Consequently, particle 2' appears in the nuclear process with the corresponding 
large wave number and the probability of the nuclear process drastically increases. -F23(^) 
has to be compared with F23{ki{Ei)) , which is the square of the Coulomb factor of the usual 
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first order process, where Ei is the kinetic energy in the relative motion of particles 2 and 3 
in the initial state of the usual, first order process. If has eV order of magnitude then 

i^23(^) » F^sihiEi)). (21) 

From this relationship it follows that although the usual first order 2 + 3^4 nuclear process 
is very strongly hindered by the repulsive Coulomb interaction due to the extremely small 
value of F23{ki{Ei)), in consequence of the appearance of the much larger quantity F23(^) 
in the rate of the second order electron assisted process the hindering effect practically 

disappears. 



IV. HEAVY PARTICLE ASSISTED NUCLEAR PROCESSES 



In electron assisted nuclear reactions heavy, charged particles are created in the deceler- 
ating process of reaction products of the electron assisted processes (e.g. of energy of 
0.0758 MeV is created in the electron assisted d + d ^ \He reaction). The energy of these 
heavy particles may be intermediately low (of about 0.01 [Mel^]) so their nuclear processes 
have to be also considered among the accountable nuclear processes. The corresponding 
graphs can be seen in FIG. 2. The graph of FIG. 2(a) depicts a nuclear capture process 
and FIG. 2(b) shows a nuclear reaction. (The processes, where the nuclear interaction is 
followed by the Coulomb interaction, can also be neglected.) Now particles 1 and 3 belong 
to system A and particle 2 which is a free, heavy particle created in an electron assisted 
nuclear process belongs to system B. According to the applied notation, particles 2 and 3 
take part in a nuclear process and particle 1 only assists it. The different processes will be 
distinguished by the type of the assisting particle and also by the type of the nuclear process. 
In our model charged, heavy particles, such as protons (p), deuterons {d) can form system 

they may be particle 2, which are supposed to move freely in a solid (e.g. in a metal). 
The particles, that take part in the processes and belong to system A are: locahzed heavy, 
charged particles (bound, localized p, d and other nuclei) as the participants of Coulomb 
scattering (particle 1) and localized heavy, charged particles (bound, localized p, d and other 
nuclei) as nuclear targets (particle 3). The problem, that there may be identical particles in 
systems A and B that are indistinguishable, is also disregarded here. 

The calculation of the transition probabihty per unit time [wf^'j of the process can 
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(b) 

FIG. 2: The graphs of heavy particle Coulomb-assisted nuclear processes. The filled dot denotes 
Coulomb-interaction and the open circle denotes nuclear (strong) interaction. The localized par- 
ticles 1 and 3 belong to system A. The free particle 2 (e.g. p or d) belongs to system B. All the 
particles are heavy and positively charged. From the point of view of the nuclear reaction particle 
2' is ingoing, particle 3 is target and particles 4, 5, 6 are products of the process. FIG. 1(a) is a 
capture process and FIG. 1(b) is a reaction. 

be performed through similar steps to those applied for the calculation of the rate of the 
electron assisted process. The main difference is that particle 1 is heavy and localized. 
The matrix element of the screened Coulomb potential modifies as 

_ 2 (27r)^ 6^2:12:2 ipu {k2i - kij - ka^) 

Here 

iPu (K) = j (x) e^^K-'^rfx, (23) 
where ipu (x) stands for the initial, localized state of particle 1. It is supposed that 

V^H(x)=(f)'''e-*-^ (24) 

is the wave function of the ground state of a 3-dimensional harmonic oscillator of angular 
frequency uji with /3i = ^JmiUJi/h 7]. The calculation of the total rates of the heavy particle 
assisted nuclear processes can be found in Appendix II. 
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The total rate of the leading, heavy particle assisted nuclear capture processes valid in 
the case of fee metals 

(«) = ^^7=Gs {k2. 2hf) hl„^,x (25) 



Here the parameter e.g. in the case of a localized deuteron, is = ^Jniduojh, where 
TUd is the deuteron rest mass, hw corresponds to the energy of an optical phonon in the 
deuterized metal, and d is the lattice parameter, 

< {a) = Ko (a) / {2nf' . (26) 

For Kq (a) see (l20l) . The quantity 

^^^^''^-exp[2vrr,,^3(A)]-l ^''^ 

with the parameter 



r,!^^(A) = Z2Zsafa23\l^^, (28) 
where mgc^ = 931.494 [MeP^] is the atomic mass unit, and 

a,fc = -^^^ (29) 

is the reduced mass number of particles j and k of mass numbers Aj and Aj.. 

It can be seen from flTTI) and f l25]) that this process opens the door for a great variety of 
nuclear processes. To get an order of magnitude estimation of the effect we take ~ A3 ^ 
A2, thus = A2 and 014 = A4/2 ~ A3/2 = ^3. With these approximations 



2Tivkapp = 2vr^2«/A2 J'l^. (30) 

We take as a typical value Aapp = 4 [MeV^] in (130|) that yields 2717723 ^^^^ = 0.497 x 2:2^2-23''^. 

If 27r?723,app — 15, one can obtain F2^3(A) > 4.59 x 10~^, that can be produced with any pair 

1/2 

of heavy particles 2 and 3 of Z2A2Z-^ < 30. From this condition one can draw an important 
and surprising conclusion: it allows a great variety of nuclear processes, that were thought 
improbable up till now. Moreover, in many cases 2717723 ^ consequently F^^{A) = 1 (see 
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(127|) ) and the hindering role of F23{A) disappears. In order to show the capabihty of the 
heavy particle assisted nuclear processes, some cases of the proton assisted proton captures 

iX + p^ iXlY + A (31) 

are investigated in Appendix III. 

V. CHARGED PARTICLE ASSISTED dd AND pd FUSION PROCESSES IN 
SOLID METALS 

As an example, let us take the simplest charged particle processes, the usual nuclear 
fusion processes: 

d + d^ ^i/e + 7 + 23.847 MeV, (32) 
d + d^ p+ + 4.033 MeV, (33) 
d + d^n+ ^He + 3.269 MeV, (34) 

p + d^ ^He + j + 5.493 MeV. (35) 

The coefficients of the astrophysical factor (see Introduction) of these nuclear fusion processes 
are available pj] . Owing to their astrophysical importance the low energy range of the above 
fusion processes has been extensively investigated. In the extremely low energy range, i.e. 
at near room temperature, the S{E) = S{0) approximation is valid and the ratio of the 
rates of the processes fl32l) . fl33|) and flMl) is determined by their 5(0) values and therefore 
it can be considered energy independent. The corresponding S'(O) values of processes (l32l) . 
[33D, dMl) and ([35]) are: 5.0 x 10-^ 5.6 x lO'^, 5.5 x lO'^ and 2.0 x 10"^ in [MeVb] units 
l|. Consequently, in a deuteron plasma or in a proton-deuteron mixed plasma the leading 
processes are the processes ( l33l) and ( IMl) with approximately the same rates, and the rates 
of processes (l32l) and (l35l) are many orders of magnitude smaller. 

For comparison we have calculated the 5'(0) values of the above processes in the simple 
nuclear model used in this article and supposing magnetic type transition in reactions (|32|) 
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and ([35]). The results are: 5.2 x 10-^ 5.0 x 10-^ 4.4 x lO'^ and 4.5 x lO'^ in [MeVb] units 
corresponding to reactions f l32|) . fl33|) . f l34|) and f l35|) . respectively, indicating that our simple 
model is able to give qualitative or semi quantitative conclusions. 

Now we review the ratios of the rates of the processes obtained in solid metals. In the 
case of the second order, simple electron assisted nuclear processes we have obtained for 
the relevant quantities (of the corresponding rates, see (HHjl . (168|) and (ITB ) -F23(^) = 0.356, 
K2K^J22 = 1-38 X 10""^ and K2J32 = 5.79 x 10""^ for the electron assisted versions of 
processes (|32|), (|33|), (IMj), respectively, and F23(^) = 0.0282 for the electron assisted version 
of process ( l35l) (see Appendix I.). From these numerical values one can conclude that, 
contrary to the above, in the family of second order electron assisted dd processes the 
electron assisted version of process ( 132|) is leading, and the electron assisted versions of 
processes (|33|) and flM|) have much lower rate. The total rate of the electron assisted version 
of dMD is 

l^i? («) = 0.012 X (£;r'^'[eV]) uN2[s-'], (36) 

— 1/2 

where Ei is the energy of the initial free electron in the conduction band. Averaging 

by means of the Fermi-Dirac distribution in the Sommerfeld free electron model at T = 

yields (^E^^^^ [c^]^ = 2 (^Ep^^'^ [^^ij? where Ep denotes the Fermi energy 8]. 

In the case of heavy particle assisted processes we have obtained F2^{A) = 0.902, and 
^^^) _ Q_QQ2g X GsN2U [s-^] for process d + d ^ ^He in the case of a deuterized Pd 
target, i.e. the particles 1, 2 and 3 are all deuterons. In the case of process p + d ^ ^He and 
in the same target material, i.e. particles 1 and 3 are deuterons, and particle 2 is a proton 
F^^{A) = 0.846 and W^^^i'^ (a) = 0.378 x GsN2U [s'^]. 

From the above rates one can conclude that if energetic, heavy charged particles are 
present in the sample the heavy particle assisted processes are not negligible and among 
all the charged particle assisted processes discussed here the electron assisted version of the 
( l32l) process is the leading one. 



VI. SUMMARY 



It is found that, contrary to the commonly accepted opinion, in a solid metal surrounding 
nuclear reactions can happen between heavy, charged particles of like (positive) charge of 
low initial energy. It is recognized, that one of the participant particles of a nuclear reaction 
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of low initial energy may pick up great momentum in a Coulomb scattering process on a 
free, third particle of the surroundings. The virtually acquired great momentum, that is 
determined by the energy of the reaction, can help to overcome the hindering Coulomb 
barrier and can highly increase the rate of the nuclear reaction even in cases when the rate 
would be otherwise negligible. It is found that the electron assisted d + d ^ ^He process 
has the leading rate. In the reactions discussed energetic charged particles are created, that 
can become (directly or after Coulomb collisions) the source of heavy charged particles of 
intermediately low (of about a few keV) energy. These heavy particles can assist nuclear 
reactions too. It is worth mentioning that the shielding of the Coulomb potential has no 
effect on the mechanisms discussed. 

n 

Our thoughts were motivated by our former theoretical findings [9] according to which 
the leading channel of the p + d ^ ^He reaction in solid environment is the so called solid 
state internal conversion process, an adapted version of ordinary internal conversion process 
[lol |. In the process formerly discussed if the reaction takes place in solid material, in 
which instead of the emission of a 7 photon, the nuclear energy is taken away by an electron 
of the environment (the metal), the Coulomb interaction induces a p + d ^ ^He nuclear 
transition. The processes discussed here can be considered as an alternative version of the 
solid state internal conversion process since it is thought that one party of the initial particles 
of the nuclear process takes part in Coulomb interaction with a charged particle of the solid 
material (e.g. of a metal). 

There may be many fields of physics where the traces of the proposed mechanism may have 
been previously appeared. It is not the aim of this work to give a systematic overview these 
fields. We only mention here two of them that are thought to be partly related or explained 
by the processes proposed. The first is the so called anomalous screening effect observed in 
low energy accelerator physics investigating astrophysical factors of nuclear reactions of low 
atomic numbers HI]. The other one is the family of low energy nuclear fusion processes. 
The physical background, discussed in the Introduction and in the first part of Section V., 
was questioned by the two decade old announcement jlj] on excess heat generation due to 
nuclear fusion reaction of deuterons at deuterized Pd cathodes during electrolysis at near 
room temperature. The paper [l^ initiated continuous experimental work whose results 



were summarized recently 



The mechanisms discussed here can explain some of the 



main problems raised in [l3|. (a) The mechanisms proposed here make low energy fusion 
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reactions and nuclear transmutations possible, (b) The processes discussed explain the lack 
of the normally expected reaction products. 

The authors are indebted to K. Hartlein for his technical assistance. 



VII. APPENDIX I. - RATE CALCULATION OF ELECTRON ASSISTED NU- 
CLEAR PROCESSES 

For particles 1 and 2 (electron and ingoing heavy particle) taking part in Coulomb in- 
teraction we use plane waves. Thus the Coulomb matrix element is calculated in the Born 
approximation 

cb 2 (27r)^ e^ziZ2 5 (ks^ + ki^ - kij - ka^) . . 

" |k,,-k,|VP ^^^^ 

which is corrected with the so called Sommerfeld factor gs (see ([9])). In the intermediate 

state we use plane waves of wave vector k2^ for particle 2. The final state of the electron is 

also a plane wave. The Dirac delta 5 {^^2^ + kij — kij — k2^) will result k2^ = k2j + ki^ — kij 

in the intermediate state (after integration over k2^, see later), i.e. in the intermediate state 

particle 2' may have large momentum, which is determined by the reaction energy due to 

the final wave vector of the electron. 

When calculating the matrix element of the strong interaction potential between particles 
2 and 3 we use the approximate form of <y9(r) given in (jG]). 

For the process 1(a) the Weisskopf approximation is used, i.e. for the final nuclear state 
of one nucleon (of particle 4) we take 



if r < i?, where R is the nucleon radius, and $ jw (i") = for r > i?. In evaluating V^* the 
long wavelength approximation (exp {i^^2^J^ ■ x) = 1) is used with sR = 1 that results 



12nR ( ^ 2 

e 



V}-'- («) = -2/^23 (^2,)^-^ ( 1 - 7 ) (39) 



in single nucleon approach. In the case of process 1(b), for the final states plane waves are 
assumed producing 

Vf;'^(/3) = -2/V23(M^ (40) 
with a neutron as one of the particles 5 and 6. If both particles 5 and 6 have positive charge 
then the expression of Vf^^ (/3) must be multiplied by /sg. The wave vectors kij, k2i of the 
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particles 1 and 2 in their initial state can be neglected in the calculation since |kij| <^ |ki/|, 
|k2i| ^ |ki/|, furthermore A ^ |ki/|. Collecting everything obtained above and applying 
~^ I [^/ (27r)^] dk2^ correspondence 



T^fitype) = B{type)z^Z26 {k^f + K^type)) x (41) 



9sf23{k2p,) 5(ki/ + k 



2m. 



{E^ - Ei) Ik 



72 '^^2^ 



1/1 



where 



B(a) = B,^il-l], (42) 



with 



So = -8(27r)'aj (^Qj {hcf. 



and 



B{nP) = i?o^, (43) 



B{+P) = B{n/3)U (44) 



For more precise result, beyond the Weisskopf and long wavelength approximations, and 
beyond the single nucleon approach the integrand of ( 14T1) must be multiplied by a model 
dependent correction factor 

Vflitype^ 

of the k — th target particle, where V^^ [type, k^) is the nuclear matrix element calculated in 
an other model, without the long wavelength approximation and beyond the single nucleon 
approach. 

We will mainly treat proton (and deuteron) capture processes in which the interaction of 
proton (deuteron) takes place with more than one nucleon. In the Weisskopf approximation 
the sum of the matrix elements of proton-proton and neutron-neutron interactions can be 
neglected due to the presence of exchange terms, so in the case of proton capture the matrix 
element of the strong interaction must be multiplied by the number of interacting neutrons 
in the nucleus. In our case it means the neutron number of particle 3 = — {hcorr,3 = 
N^). In the case of deuteron capture reactions and also in the Weisskopf approximation the 
matrix-element must be multiplied by ^3 (/icorr.s = ^3)- 
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A. Rates of electron assisted p or d capture processes 

First the case type = a is treated. Neglecting the initial kinetic energies and £^2^* in 
the denominator of flT^ . — Ei = Eif + E2^ Eif that will have a value Eif A 
because of the energy Dirac delta. Using the ~^ I [^/ (27r)^] c?k2^ and the J2f ~^ 
J [V/ {2TTf]d'ki X J [V/ {2nf]d'kif substitutions, the identities [{2nf 6 {k^ + kif)]^ = 
(27r)^5(k4 + ki/)(27r)^5(0) and {2nf 5 (0) = V, and the hfdhf = EifdE^f/ {hcf rela- 
tions, furthermore carrying out integrations one can obtain 

(a) = Ko (a) Gs {ku, 2hf) ^f^h^^^^,. (46) 

Here Gs {kii,2kif) = [see ( ITOj) ]. remember that Fj^ik) = ffi^ik) [see ([8])] and Kq^o) is 
defined in ( !20l) . In the case Eij ^ nieC^ 

kif = A/ {he) . (47) 

If one is interested in the total rate W^J (a) of a sample then Wj- must be multiplied 
by the numbers of initial particles, A^^i, N2 and corresponding to electrons, and particles 
2 and 3, respectively. The quantity Ni/V = Qe/vc, where Vc is the volume of the elementary 
cell of the solid metal and Qe is the number of the valence electron states corresponding to 
one unit cell, e.g. ge = 10 and Vc = d^/i in the cases of Ni and Pd will be discussed later. 
We introduce N^/V = n-^, the number density of particles 3. It is reasonable to take for the 
number density of the target = 2u/vc = u8/d^ in the case of fee metals (such as Pd 
and Ni) with u denoting the deuteron (or proton) over metal number densities. The result 
is ^ 

W}^^ (a) = K,,t (a) {Gsiku, 2k,j))^^ E^^hl^^^^uN, (48) 

where Ktot («) is determined by (fT9|) and the average is made by means of Fermi-Dirac 
distribution. 

B. Rates of electron assisted processes with two outgoing heavy particles 

In the case of type = +/3 or nfi (for the process (6) of FIG.l) the calculation is modified 
as follows. Let the wave vectors of particles 5 and 6 be ks and kg. This time, the — >■ 
/ \V/ (27r)^] rfks X / \V/ i2Tif] rfke x / \V/ (27r)^] dki/ correspondence and the identities 
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[5 (ki/ + kg + kg)]^ = 5 (ki/ + kg + kg) 5 (0) and (27r)^ S (0) = V are used. The integration 
over kif results kij = — ks — kg. Using again the N^/V = = 2u/vc = u8/cP relation, 

Wiol {(3) = Ktot (a) 2N2 ( ^ J ^^corr,3W^2 (49) 



with 



j-i L , , , x (50) 

J y |/t5 + Kel e{j + Ke) 



where 



S{e) = 5 [es (k^) + ee (kq) + eif (^5 + kq) -Se-l]. (51) 

Here kj = /ickj/A and €j = Ej/A = (2mj)] A are dimensionless quantities with Ej as 
the kinetic energy of particle j (j = 5 or 6), 5e = rneC^/ A, Sif = Eif / A = \J + kq\^ + ^l — 
5e and the suffixes 23 of F23 and 56 of F56 refer to particles 2, 3 and 5, 6, respectively. If 
one of the particles 5 and 6 is a neutron then 

i^56 = 1. (52) 

It is useful to introduce the following new variables in fl50l) : 

a = + (53) 

and 

b = - Ke. (54) 
Thus, if both particles 5 and 6 are positively charged then J = J^j^, 

J+p = A23A5GJ2 (55) 

where 

A23 = 27iZ2Z3affi23c'^/A, (56) 
= 2Trz5ZQaf^5ec^/A (57) 

and 

f A23\ ( A56\ , „ 

X exp I — — - exp — -rr-j- aactb. 
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If one of the particles 5 and 6 is a neutron then J = Jnp, 

JnP = A23J3 (59) 

with 

Here and above 

6 [e (a, b)] = 6 ( ^^b' + ./^^Tsl - 1 - 4 V (61) 

In obtaining (ICTj) the Aa^/ (S/isgC^) and A (^ ^^^^2 — ^mec'^ ) ^ ' ^ terms in the argument can 
be neglected since A/ (/isec^) ^ 1, ct ^ 1, and the dominant range in the integrals is where 
6 ^ a. Using the 6 [g (x)] = 5 (x — Xi) / g' (xi) identity, where Xi is the root of the equation 
g {x) = 0, the integrals J2 and J3 can be written as 



with 



J2 = P2J22 (62) 



132 = 647r^/i56c7A, (63) 



exp f— ^) exp ( — ^ 



A56a/A/(8/JS6C^) 



and 



with 



J3 = /33 J32 (65) 

/33 = 5127r2(^56cVA)', (66) 

/•an^ax exp (-^) f 1 + 4 - v^a^T^f) rfa 
■^32 = / "—^ ,2 ^ (67) 

where Omax = + Using all the above results, the total rate W^J (+/3) of the process 
having two charged, heavy products reads 

Wi^} (+/3) = Ku., (a) K2K^ {Gs (h,, 2hj))^^ :^hl„^,uN2, (68) 

where 

512 I R\ 21 2\2 lnr^\ 

^ ^^^^ 3(i_ 2^2 ^ \nc ) "^^^^^ ^^^^^ ^ ^ 
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and 



TT 



(70) 



(2') 

The total rate Wi^l {nf3) of the process, in which one of the two heavy products is a neutron, 
reads 



J: 



32,2 



(71) 



In order to compare the total rates psj) . (!68|) and (!7T|) of the processes of different type, 
the quantities -F23(^), K2K+J22 and K2J32 have to be compared. The effect of factor 
{Gs {kii,2kif)) is not essential in the electron assisted processes. 

VIII. APPENDIX II. - RATE CALCULATION OF HEAVY PARTICLE ASSISTED 
NUCLEAR PROCESSES 



In the case of heavy particle assisted nuclear processes the quantity Tj^^^ modifies as 



Tlf{type) = B{type)ziZ26 (ki/ + K(^type)) 



X 



(72) 



gsf23{k2fi) ^li (k2i - ki/ - ka^; 



{Ef^ — Ei) |k 



where 



E, 



2i ~ '^2fi\ 



,2 ,2 

""^1 ,/ + ~ ko 



+ A2 



c/k 



2/i 



EM = B, 



with 



5(n/3) 



(l-f)V 31/ 



B{a) UttR^ 



and 



Since /3i < IKI 



B{+f3) = B{nP)U 
iif\, the Fourier transform of 

V^i.(K) 



23/2^3/4 _K 



(3t 



3/2 



-e ^"1 



(73) 

(74) 

(75) 
(76) 

(77) 
(78) 
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allows the approximation 

MK) = 87r9/^/3f '5(K) (79) 



in f l72|) . As a result the integral over in fl72|) can be carried out, while |k2j| and A 
can be neglected since |k2j| ^ |kij| and A ^ |ki/|. Neglecting Ei in the denominator, 
E^ — Ei = fi^klj/ {2^12). The result for Tlj\type) obtained in this way is 



Tlf\type) = B{type)ziZ26 {kif + K^typ,)) x (80) 

lb A_f 1 f 



A. Rates of heavy particle assisted nuclear capture processes 



First the process of type = a is dealt with. Substituting ( I8OII into ( ITTj) the 

/ 



J2f ^ I \yl (27r)^] rfk4 X / \yl (27r)^] dki/ correspondence, the identities [b (ki/ + k4)]^ 



b (ki/ + k4) b (0) and (27r)'^ 5 (0) = ^ are used. When integrating over k4, the b (ki/ + k4) 
leads to the k4 = — ki/ replacement. Therefore Ef = Eif (/ci/) + -£^4 (A^i/) = ^"^kfj^/ (2/ii4) in 
the energy Dirac delta. Applying the kifdkif = ^udEf/h"^ and kif = y^2ji^^Ef /h relations 
the integral over kif is converted to an integral over Ef and it is carried out with the aid 
of the energy Dirac delta. Thus the transition rate of a heavy particle assisted reaction of 
type = a 



'ai2y {z^Z2fF^%{A) 1 



^ 2),. ^ flo^Gs {k2., 2k,f) hl„,x (81) 



X 



aij y/aU V' 

For -ft^Q (a) see (!26l) . the quantity ^2^3 (A) is given by (!27j) . the parameter 7733 (A) is determined 
by f l28|) . moc^ = 931.494 [Me\^] is the atomic mass unit, and ajk is the reduced mass number 
of particles j and k of mass numbers Aj and given by ( l29l) . 



The total rate w/^^^^''' (a) can be obtained multiplying wjf''^ (a) by A2 and A''i3. A^2 is the 
number of particles 2 and Ni^ is the number of possible pairs of particles 1 and 3. Introducing 
ni3 = N13/V, the number density of particles 1 and 3, we take again rirs = 2u/vc = 8u/d^ 
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valid in the case of fee metals. 



(«) = (fc,., 2k,j) hl„,x (82) 



auj JoYa A^/^ 



The parameter /3i in the case of a localized deuteron is /3f = radio jfi^ where is 
the deuteron rest mass. In Pd, that is one of the possible materials which can absorb 



hydrogen isotopes well, the energy of the ground state of the oscillator Eq = ^hu = 72 
[meV] jl^ leading to hco = 48 [meV], that corresponds to the energy of an optical phonon, 
and /3f = 4.81 x 10^ [cm^^]. This value alters as /3f = 4.04 x 10*^ [cm"^] for protons. With 
these numbers the characteristic constants, which appear in W^^^'^^ (a) ( l82l) . are 

and 

d^y/niQC^ 

in the case of deuterized and protonated Pd, and in the cases when particle 1 is a deuteron 
or a proton, respectively. 



190 [MeV'^/h-^] (83) 



113 [MeV'^/h-^] (84) 



B. Rates of heavy particle assisted processes with two outgoing heavy particles 

In the case of type = +/3 or n[5 (for the process (b) of FIG. 2) the calculation is modified 
as follows. Now the E/ ^ / [^1 (27r)^] c/kgx / [V/ {2Txf] c/kex / [V/ (27r)^] rfki/ correspon- 
dence and the identities [5 (kij + ks + kg)]^ = 5 (kij + ks + kg) 5 (0) and (27r)^ 5 (0) = V 
are used and the integration over kij leads to the kij = — ks — kg replacement. In the 
argument of the Dirac delta = (jks + kg]) +-E5 (^5) +-^6 (^e); which can be rewritten 
as 

Ef = -^ + -^ + —hkeC, 85 
2/ii5 2/ii6 mi 

where ( = cosGsg with 656 the angle of vectors ks and kg. Now new variables 



Xj 



fi2 

k„ J =5,6 (86) 



2miA 



are introduced. With these variables 



Ef = { —xj + —xl + 2x5XeC ) A (87) 
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and 

^1/ = H + ^6 + S^s^^eC) ■ (88) 

Applying the 6 [x (C)] = \x' (Ci)!"^ ^ {C ~ Ci) identity, where is the root of equation x (C) = 
0, the integration over ^ can be carried out. In our case 

ci = ^ ■ (89) 

Introducing the notation 

g[±) (2^5, xq) = + ^6 ± 2X5X6C1 (90) 

and after some algebra the transition rate 

<^'' ^^yp^) - ^r^Gs ih., 2h,) /.L.,3X (91) 

\2 



X {Z1Z2) { ) JtypeTT: 



where 



2\ 2 



Xo(^)=2V/24(^Q i?%c^. (92) 
For type = +/3 the integral Jtype is 

J+/3 = / / ^(63)^(^56) x^dx^xedxQ (93) 
Jo Jo i/+(a:^5,a:;6) 

and for type = the integral J^j^pe is 

Jnp= / ^7 -xscfxgxedxe. (94) 

9+{X5,X6) 



Jo Jo 

Here the following notation is used: 



and the variables are 



and 



63 = '^^^^^^^f^^=^9(+)^ {xb, xe) (96) 



with 



g{±) {X5, xe)^xl(lT—]+xl(lT—)±l. 



(98) 
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The upper limits x^^up and x^^up are determined by the condition < 1. 

Now the total rates wj;^]''^ {type) of reactions type = +/3 or njH can also be obtained 

if) 



multiplying wif''^ (type) by N2 and A''i3, and using nis = Ni^/V = 8u/d^ 



W,T {type) = ^^^j^Gs ih,, 2ky) hi^,^,x (99) 



2 fP'12C^\ 
X {ziZ2f ( ) JtypeUN2. 

In the rates W^'^}'^ [a) ([H2D and W^'^}''" {type) ^ the function Gs{k2i,2kif) = 
Fi2{k2i) / Fi2{2kif) plays important role since Fi2{k2i{Ei)) has strong energy dependence. 

— 1/2 

If e.g. both particles 1 and 2 are protons then 27r?7i2 {E2i) = 0.700 x £2^ resulting 
Fi2{k2;{E2i)) = 6.4 X 10-3at E2i = 0.01 [MeV]. 

It is also an important aspect of the mechanism that if particle 1 is a proton then in 
each proton assisted nuclear capture process an energetic proton (of energy of about a few 
MeV) is created too. This proton in its decelerating processes can create (secondary) free 
protons in the crystal, if localized protons are also present. The secondary free protons 
may take part in further proton assisted nuclear capture processes. Thus in this process the 
secondary protons can play a role similar to that played by the secondary neutrons in the 
case of nuclear fission. 

IX. APPENDIX III. - PROTON ASSISTED PROTON CAPTURE PROCESSES 

Here some cases of the proton assisted proton captures 

^X + p^ ^XlY + A (100) 

are investigated. Particles 1 and 2 are protons. Particles 3 and 4 have mass numbers A and 
y4 + 1, respectively (see FIG. 2). As a first example we take protonated Ni. In this case the 
possible processes are 

^m + p^ ^+^Cu + A. (101) 

Most of the daughter nuclei 29'^Cu decay by the 

i+^Cu + e^ i^^m + QEc (102) 
electron capture reaction. TABLE I. contains the relevant data for reactions ( llOip and ( I102p . 
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A 58 59 60 61 62 64 

A(MeV) 3.417 4.479 4.801 5.867 6.122 7.452 

TA 0.68077 (2.4 x lO^^s) 0.26223 0.0114 0.03634 0.00926 

27rr?^3 15.04 13.14 12.69 11.48 11.24 10.18 

sa(A) X 10^(MeV-'^/2) 4 _ 4 23 0.276 0.945 0.316 

t{s) 81.5 1422 1.20 x 10^ 584.4 stable stable 

Q£;c(MeV) 4.789 6.127 2.237 3.948 



TABLE I: Numerical data of the ^^Ni + p ^ 2^^Cu + A and 29 ^Cu + e ^ + Qec 

reactions. Computed values of 2TTr]2^ and sa (A) are given (for their formula see the text). A is 
the mass number, is the relative natural abundance, r is the half life of the 29^^Cu isotope 
produced. It decays by electron capture of reaction energy Qec- In the row ta 3X A = 59 the 
2.4 X 10^^ s stands for the life time of the unstable Is-^^- 

An other interesting family of the heavy particle (proton) assisted proton capture is 

i^Pd + p-^ ^7+1^(7 + A, (103) 

that is mainly followed by the 

^+iyl^ + e^ itPd + QEC (104) 
reaction. (The relevant data can be found in TABLE II. The nuclear data of TABLES I. 



and II. are taken from 15|. 



For the processes fllOip and f llOSp ai2 = 1/2 and an = 023 = 1 is a good approximation 
in the case of protonated Ni and Pd. First we take a protonated Ni {z^ = 28). The rate 
^tot''^ (a, A) belonging to process (11 01 1) is 

WS'' («, A) = G.^^^^M4£ {A - z,f SA (A) uN2. (105) 

Here Kq (a) is given by ( 126|) . d is the lattice parameter, moc^ = 931.494 [MeK] is 
the atomic mass unit, and for /3f see Appendix II. A under (|82|) . These all result 
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A 


102 


104 


105 


106 


108 


110 


A(MeV) 


4.155 


4.966 


5.814 


5.789 


6.487 


7.156 




0.0102 


0.1114 


0.2233 


0.2733 


0.2646 


0.1172 


27r77t3 


22.41 


20.49 


18,94 


18.98 


17.93 


17.07 


SA (A) X 10"(MeV-^/2) 


0.0285 


1.06 


5.31 


6.35 


11.1 


7.80 


t{s) 


4002 


3.56 X 106 


1438 


stable 


stable 


6.44 X 10^ 


Qi,c(MeV) 


2.688 


1.346 


2.965 






no EC 



TABLE II: Numerical data of the ^^Pd + p fr'^Ag + A and ij^'^Ag + e t^^Pd + Qec 
reactions. Computed values of 27rr/23 and sa (A) are given (for their formula see the text). A is 
the mass number, is the relative natural abundance, r is the half life of the f^^^Ag isotope 
produced. It decays by electron capture of reaction energy Qec- 

2K^ {a) {0[f d-^ (moc^)"^/^ = 28.3 [Mel^^/^^-ij _ jy^ jg ^he actual number of free protons 
in the sample, 

.^(A)=r^^, (106) 

where is the relative natural abundance and 

, 27.8 [Mel/V^i 
2^4, = L_ 1 (107) 



with zs = 28 in (!28|) . The formula f llOSp can be also used in the case of processes (11031) with 
Sa and A data of TABLE II. and using 



^ _ 45.7 [Me1/i/2] 



= (108) 



that comes from with = 46 (corresponding to Pd). 

The processes (llOip and (I103P may produce stable 29CU, %Cu and "^^^Ag, 47^^(7 isotopes, 
respectively, whose heavy particle assisted proton capture reaction may give rise to a chain 
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of nuclear transmutations. 



[1] C. Angulo et al, Nucl.Phys. A 656, 3-183 (1999). 

[2] P. Marmier and E. Sheldon, Physics of Nuclei and Particles, Vol. 1 (Academic, New York, 

1989) pp. 315-317. 
[3] K. Alder et al, Rev. Mod. Phys. 28, 432-542 (1956). 

[4] W. Heitler, The Quantum Theory of Radiation, 3rd ed. (Clarendon, Oxford, 1954) Ch. V.,25§, 
(19). 

[5] W. H. Louisell, Quantum Statistical Properties of Radiation (Wiley, New York, 1973) pp. 
81-85. 

[6] J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics (McGraw-Hill, New York, 
1964). 

[7] C. Cohen-Tannoudji, B. Diu, and F. Laloe, Quantum Mechanics, Vol. 1. (Wiley, New 

York/English version/, Hermann, Paris, 1977). 
[8] J. Solyom, Fundamentals of the Physics of Solids, Vol. H., Electronic pro;?eriies (Springer, 

Berlin-Heidelberg, 2009), pp. 35-37. 
[9] P. Kalman and T. Kesztlielyi, Pliys. Rev. C69, 031606(R) (2004); Phys. Rev. C79, 031602(R) 

(2009). 

[10] J. H. Hamilton, Internal Conversion Processes (Academic, New York, 1966). 

[11] F. Raiola et al, Eur. Phys. J. A 13, 377-382 (2002); Phys. Lett. B 547, 193-199 (2002); 
C. Bonomo et al, Nucl. Phys. A719, 37c-42c (2003); J. Kasagi et al, J. Phys. Soc. Japan, 
71, 2881-2885 (2002); K. Czerski et al, Europhys. Lett. 54, 449-455 (2001); Nucl. Instr. and 
Meth. B 193, 183-187 (2002); A. Huke, K. Czerski and P. Heide, Nucl. Phys. A719, 279c-282c 
(2003); A. Huke et al, Phys. Rev. C 78, 015803 (2008). 

[12] M. Fleishmann and S. Pons, J. Electroanal. Chem. 261, 301-308 (1989). 

[13] E. Storms, Naturwissenschaftcn, 97, 861-881 (2010). 

[14] E. Wicke, H. Brodowsky, and H. Ziichner, Hydrogen in Palladium and Palladium Alloys, in 
Hydrogen in Metals, edited by G. Alefeld and J. Volkl, Vol. 2., Ch. 3. (Springer, Heidelberg, 
1978) Fig. 3.7 p. 88. 

[15] R. B. Firestone and V.S. Shirly, Tables of Isotopes, 8th ed. (Wiley, New York, 1996). 

27 



